Bottonium mass - evaluation using renormalon cancellation 
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We present a method of calculating the bottonium mass Mr (IS) = 2m b +E bb . The binding energy is separated 
into the soft and ultrasoft components E bb = E bb (s) + E bb (us) by requiring the reproduction of the correct residue 
parameter value of the renormalon singularity for the renormalon cancellation in the sum 2m b +E bb (s). The Borel 
resummation is then performed separately for 2m b and E bb (s), using the infrared safe m b mass as input. E bb (us) 
is estimated. Comparing the result with the measured value of Mt(1S), the extracted value of the quark mass 
is rn b (jj, — fn b ) = 4.241 ± 0.068 GeV (for the central value a a (Mz) = 0.1180). This value of m b is close to the 
earlier values obtained from the QCD spectral sum rules, but lower than from pQCD evaluations without the 
renormalon structure for heavy quarkonia. 



Heavy quarkonia qq (q — 6, t) can be inves- 
tigated by perturbative methods (pQCD) via ef- 
fective theories NPQCD Q] and pNRQCD g] (or: 
vNRQCD |3]) because of the scale hierarchies of 
the problem: m q > m q a s (/x s ) > m q a 2 (iJ, us ) ~ 
Aqcd- Here, m q is the (pole) mass of the quark, 
fi s ^m q a s (^ s ) is the soft, and fi us ~ m q a 2 (fj, us ) 
the ultrasoft energy. The quarkonium mass is 
M q q = 2m q + E qqi where the binding energy 
consists of the soft and ultrasoft regime contri- 
butions: E qq — E qq (s) + E qq (us). A practical 
problem which appears in the course of evalua- 
tion of M qq is that the perturbative pole mass has 
an inherent ambiguity Sm q ~ Aqcd (~0.1 GeV) 
due to the infrared (IR) renormalon singularity 
which appears at the value of the Borel transform 
variable b — 1/2 for m q /m q . Here, rn q is the in- 
frared safe (renormalon-free) MS mass. However, 
the static potential V qq (r) has a related ambiguity 
SV qq (r) - Aqcd such that 8(2m q + V qq ) = 0, i.e., 
the renormalon singularity cancels for the com- 
bined quantity 2m q + V qq 0] (see also |H]). The 
static potential is a quantity which does not con- 
tain ultrasoft regime contributions [Hj- Therefore, 
Eqq(s) contains the entire V qq and kinetic energy 
effects, the latter are renormalon-free. Thus, the 
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6=1/2 renormalon singularity of V qq and E qq (s) 
are equal and hence the singularity must cancel 
also in the combination 2m q + E q q(s) 

S[2m q + E m (s)} =0 . (1) 

In principle, E qq (us) could be included in this re- 
lation. However, in practice, it distorts the can- 
cellation effects since we know these quantities 
only to a hnite order in perturbation expansions 
[cf. the discussion following Eq. {T2|l]. 

One possibility of evaluating the bottonium 
ground state T(1S) mass is to use an infrared-safe 
(renormalon-free) quark mass (mb, m^ s , etc.) 
and a common couplant a(/x) = a s (/j,)/ir as in- 
puts in the evaluation of the available truncated 
perturbation expansion (TPS) for Mr (IS*) = 
2mb + E b i, in order to avoid the 6=1/2 renor- 
malon (divergence) problems throughout [7], and 
then extract the value of mb from the measured 
value M T (IS) = 9460 MeV. 

Another possibility is to evaluate E b b in terms 
of the pole mass m& and of a (/it), taking into 
account the 6 = 1/2 singularity of E bb (using, 
e.g., the Principal Value [PV] prescription in the 
Borel integration), and adding 2uib to E bb . From 
M r (lS) = 2mj, PV) +E bb (m { b PV) ), the PV-value 
of the pole mass m b is then extracted, and subse- 
quently the value of rn b (via PV Borel integration 
prescription) . This is the approach of Ref. [Sj . 

Our approach _9 follows to a significant de- 
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gree the latter approach, but with some impor- 
tant modifications: 

1. The input parameter is the renormalon-free 
mass nib = = m b ) (and, of course, 
the QCD couplant a(/i)). The pole mass 
m b = mftfrnj; ji m ) is evaluated via Borel in- 
tegration, accounting for the 6=1/2 sin- 
gularity, and using a hard renormalization 
scale /i m ~rnfe. The residue parameter N m 
of the 6 = 1/2 singularity is evaluated from 
the available TPS for mb/rnt,. 

2. On the basis of the knowledge of N m , we 
separate the binding energy into the soft 
(s) and ultrasoft (us) regime contributions: 
E bb = E bh( s ; Vf)+ E bb( us ; Vf), where the s- 
us factorization scale (if parametrizes the 
separation. The separation is performed 
by accounting for the renormalon cancella- 
tion in the sum 2m b + E bb (s): N m (m b ) = 
N m (E bb (s] fJ>f))- The latter relation fixes \if 
and thus the separation. 

3. The soft binding energy E bb (s\\if) is then 
evaluated via Borel integration, accounting 
for the 6=1/2 singularity (using the same 
prescription, e.g. PV, as for m b ), and using 
a soft renormalization scale /i s ~mf,a s . 

4. The value of the ultrasoft part E bb (us; fif) 
is estimated. 

5. From 2m b + E bb (s) + E bb (us) = My (IS), 
the value of m b is extracted. 

For details, we refer to Ref. 

1. Evaluation of m b and N m 

This part has been performed mostly in 
Refs. |7I8I10| . The pole mass is known to NLO: 
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) } a J (ii m )rj(ii m ) , (2) 
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where r% and ri are known coefficients (tq = 1), 
e.g., in the MS scheme, and they depend on 
(/i m A>T&); /Lt m ~TOh. The Borel transform is 



B s (b) 



ri 



r-2 



0(b 3 



(3) 



m b (l-2b)^f^ o 



c k (l-2b) k + B^(b), (4) 



where ft = (ll-2n//3)/4, ft = (102-38n//3)/16, 
v = ft/(2^o) ( n f — 4); Co = 1 and the next three 
coefficients Ck are known ([7] for k= 1,2; |§] for 

fc = 3). fig (6) is the analytic part in the bilocal 
expansion J2J jS| , and it is known up to ~ 6 2 . The 
residue parameter N m in Eq. can be obtained 
with high precision |7I8I10| 
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where, according to J2| 

i? s (6; Mm ) = (1 - 26) 1+ ^ s (6; Mm ) . (6) 

Applying the Padc P[l/1] to the known NNLO 
TPS of R s (b) then gives 



N m (n f =A) = 0.555 ±0.020 



(7) 



The pole mass m b , with m b and a(fjL m ) as input, 
can now be evaluated by Borel integration using 
the bilocal expression |@} 



S(b) = \ Re / d6exp 



fta( / u T1 



B s (6;/i m ),(8) 



where the integration path can be taken along a 
ray in the first or fourth quadrant (the generalized 
PV prescription |11I12I13| 1. 

2. Separation 



The TPS of the binding energy E, 
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is known to the impressive order 0(m b a 5 ) |14l 
11511 fill 711 811 UET)|. i.e., in Eq. © /*, (fc = 1,2,3) 
are known (/q = 1). The renormalization scale 
used in expansion JjJJ should be soft (/j s ~m b a s ) 
or lower. The ultrasoft contributions appear for 
the first time at ~m b a 5 19 20 , i.e., / 3 = / 3 (s) + 
fs(us). The us coefficient can be written as 9 : 
f 3 (us)/ir 3 = 27.5 + 7.1 lna s (^ s )- 14.2 InK, where 
k ~ 1 is the parameter of the s-us factorization 
scale /i/: /i/ = Kmi,« s (/i s ) 3 ' 2 . It can be fixed by 
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the requirement of the renormalon cancellation in 
2m b +E bb (s): 



2tt m b a(ns) 



RF(s){b; fi a ;fif)\ 



6=1/2 



(10) 



9 n s 

where, in analogy with Rs of © 

R F(s) (b; mm) = (\~2b) 1+v B F[s) {b- p.; fi f ) , (11) 

and i?_F( s ) is the Borel transform of the quantity 
F(s) = -{§/{Air 2 ))E bb -{s)/{m b a{n. s )) [in analogy 
with 5 of ©]. Since now the TPS of R F{s) is 
known to ~ b 3 , the Pade P[2/l](6) thereof can 
be taken; using then the value Q of N m , the 
renormalon cancellation condition (|10|) gives nu- 
merically the s-us separation parameter 



k = 0.59 ±0.19 . 



(12) 



It was possible to obtain the value of [if k) 
because the dependence on [if in N m of Eq. (jTU|l 
was taken (and is known) only to the leading or- 
der, the ultrasoft part was excluded, and N m is 
well-known (0 . This is similar to the scale- fixing 
in the effective charge (ECH) method [2U • If the 
ultrasoft contributions are included in Eq. lfTU|) . 
the value of N m cannot be reproduced. 

3. Evaluation of the soft contributions 

Knowing now the expansion of F(s) = 
-(9/(ATr 2 ))E b - b (s)/(m b a(n s )) up to ~ a 4 , the 
Borel transform of this quantity can be con- 
structed, e.g., with the approach of the "<j- 
regularized" bilocal expansion [H] , which is a gen- 
eralization of the bilocal expansion (J2J 



B F{s) (b) = — 



27rm b a(^ s )(l-26) 1 - 



J2C k (l-2b) 1 
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+ 4? S )W- (13) 



The exponential was introduced in order to sup- 
press the renormalon part away from b w 1/2. 
The first four coefficients Ck are known (Co = 1), 
and the analytic part is known now up to ~ 6 3 . 
The analytic part we can evaluate either as TPS 
or as Pade P [2/1] (b) . The requirement of the ab- 
sence of the pole around b = 1/2 in that part, 
and the independence (weak dependence) on the 



renormalization scale (i s for the Borel-resummed 
result E bb (s), lead us to fix the a parameter to 
the values a = 0.36 ± 0.03. The Borel integration 
is performed as in Eq. ©, with the ray (PV) path 
prescription taken. 

4. Estimate of the ultrasoft contribution 

The ultrasoft part of the energy is known only 
to the leading order (^m b a°) 



E bb{ us ) (P) ~ -~gm q TT 2 f 3 (us)a 5 (n us ) 
« (-150 ±100) MeV . 



(14) 



Here, f 3 (us;[if) was determined in Sec. |2 the 
ultrasoft renormalization scale fi us should be ~ 
a 2 m b , but was taken numerically to be higher, in 
the soft regime (/x w 1.5-2.0 GeV =^ ct s (fi) « 
0.30-0.35), because perturbative QCD does not 
allow a running to very low scales. The bottom 
mass value was taken m b = 4.2 GeV. The non- 
perturbative contribution comes primarily from 
the gluonic condensate and gives E bb (us)( np ' « 
50 ± 35 MeV if the gluon condensate values 
((a s /vr)G 2 } = 0.009 ±0.007 GeV 4 [22] are taken. 
This then results in the following estimate of the 
ultrasoft contributions to the binding energy 



E bb (us) {p+np) w (-100 ± 106) MeV 



(15) 



In addition, there are contributions to the T(15) 
mass due to the nonzero mass of the charm quark 
123 SM r (lS, m c / 0) w 25 ± 10 MeV. 

5. Extraction of the mass rh~ b 

Adding together the Borel-resummed values 
2m b , E bb (s) and E bb {us), requiring the reproduc- 
tion of the measured mass value My {IS) (with 
the mentioned m c effect subtracted), we 
extract the following value for the mass m b = 
m b (p = m b ): 



m b (m b ) = 4.241 ± 0.068 GeV 



(16) 



when the QCD coupling value is taken as 
a s (M z ) = 0.1180 ± 0.0015. The major source of 
uncertainty in the result 1)16(1 is the uncertainty 
from the ultrasoft contributions (|15fl (±0.049 
GeV). The other appreciable uncertainties are 
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from the ambiguity of the soft renormalization 
scale fi B = 3±1 GeV (±0.013 GeV) and of 
a s {M z ) = 0.1180 ±0.0015 (±0.013 GeV). If the 
central value of the gluon condensate ((a s /Ti)G 2 ) 
is increased from 0.009 gl] to 0.024 GeV 4 (used 
in |24I25| L the central value l|16|l decreases to 
mb(jnb) — 4.204 GeV. This is close to the values 
of QCD spectral sum rule calculations which gave 
central values m b {m b ) = 4.20 GeV gj; 4.24 GeV 
[25]; and m b (m b ) = 4.23 GeV [H| (Ref. [S3 uses 
central condensate value ((a s /7r)G 2 ) = 0.019). 
The TPS evaluation of Mr (IS), without account- 
ing for the renormalon problem, extracts higher 
central values m b (m b ) — 4.349 |2"T)] . 
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